Optical response of finite-length carbon nanotubes 
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Optical response of finite-length metallic carbon nanotubes is calculated including effects of in- 
duced edge charges in a self-consistent manner. The results show that the main resonance corre- 
sponding to excitation of the fundamental plasmon mode with wave vector -n /I with / being the tube 
length is quite robust and unaffected. This arises because the strong electric field associated with 
edge charges is screened and decays rapidly inside the nanotube. For higher-frequency resonances, 
the field starts to be mixed and tends to shift resonances to higher frequencies. 



I. INTRODUCTION 

Carbon nanotubes are a one-dimensional conductor 
consisting of rolled graphite sheets^ Usually a single- wall 
carbon nanotube has a diameter of the order of a few nm 
and a length of the order of several jim. Recently, cut- 
ting and length selection of single- wall carbon nanotubes 
became possiblej^i^i^i^i^iiiSiiSiiiS In short nanotubes with 
good conductivity, edge effects start to play important 
roles in their electronic properties. The purpose of this 
paper is to theoretically study low-frequency optical re- 
sponse in metallic carbon nanotubes with finite length. 

In optical response, induced charges associated with 
optical transitions often play important roles. For exam- 
ple, effects of an electric field induced by the polarization 
of nanotubes are quite significant and tend to suppress 
interband transitions for light with electric-field polariza- 
tion perpendicular to the axis."'^^'"^^ However, strong exci- 
ton effects due to the one-dimensional featureHiM cause 
reappearance of exciton peaks even in the presence of the 
strong depolarization eSect^^ 

For electric field parallel to the axis, this depolariza- 
tion effect does not play a role and the absorption is 
described by the usual dynamical conductivity. The dy- 
namical conductivity in metallic nanotubes in the low- 
frequency region was studie d^^i^^ in relation to the ab- 
sence of backscatteringi2i2£iiSi and the presence of a per- 
fect channel.— In optical response of finite-length car- 
bon nanotubes, often called antenna effects, accumulated 
charges at both ends of nanotubes can become important 
even in field parallel to the axis. 

There have been several theoretical calculations on 
light scattering by finite-length nanotubes in low- 
frequency region. Direct analysis based on the integral 
equations of electromagnetics was reported.— i^ii^ Nu- 
merical studies were reported on electromagnetic waves 
in a planar nanotube array .^^i^^ An equivalent circuit 
model with quantum capacitance and kinetic inductance 
was considered. 

Recently, nanoscale antenna operation of a carbon nan- 
otube array was experimentally demonstrated.^^ Optical 
absorption in low-frequency (~THz) region was reported 
for nanotube bundles, where a sharp absorption peak 
was observed in aligned nanotubes but not in sprayed 



samples.— A broad absorption peak was also reported, 
but attributed to a narrow gap in quasi-metallic or 
narrow-gap nanotubes. ^-i^i^ Finite-length effects have 
been considered in connection with various phenomena 
such as Raman spectroscopy,—!^ tt plasmon absorpt- 
ion,= etc. 

In this paper, we shall calculate optical response of 
finite-length metallic nanotubes in much lower frequency 
region than an inter-band transition. In §2, the method 
to calculate the optical response is described and an 
approximation based on excitation of a single plasmon 
mode in an infinitely long nanotube is introduced. Nu- 
merical results are presented in §3. A discussion and 
summary are given in §4. 



II. FORMULATION 
A. Optical Response 

We consider a carbon nanotube with a finite length ^, 
lying along the y direction in the range —//2 < y < +1/2. 
Let i?oxt(y)e~'"* be an external electric field of incident 
light and E{y, a;)e~'"* be the effective electric field in- 
cluding effects of polarization charges. The response of 
the system can generally be described by a nonlocal con- 
ductivity (j{y,y'). 

Then, the induced current j(y)e~"^* is given by 



jiy)= o{y,y')E{y')dy'. 



(1) 



The corresponding induced charge p{y)e is deter- 
mined by the equation of continuity: 



dt dy 



= 0, 



(2) 



P{V) = -^.liv) = -|- / dy' a{y,y')E{y'). (3) 
iLu oy luj oy J 

The corresponding scalar potential (/)(?/)e~"^* becomes 

<^(y) = / K{y-y')p{y')Ay', (4) 
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cylindrical charge distributiori j^'^i'^'^'^^i'^^i^^i'^-'^i'^^ given by 



where K{y) is the kernel of the Coulomb interaction for 

ionf 

with 



2tt 



(5) 



Kiq) = -Io{\q\R)Ko{\q\R), 

K 



(6) 



where R is the diameter of the nanotube, n is the static 
dielectric constant of the environment, and In(t) and 
Kn{t) are the modified Bessel function of the first and 
second kind, respectively. The total electric field becomes 

= £;ext(y) + ^ Jdy' Jdy"K\y-y'My',y")E{y"), 

(7) 

with 

^'^y-y'^^^'^'^y-y'^^ J^K\q)e^''^y-y'\ (s) 

where 

K'iq) - Kiq)q\ 
The power absorption is given by 



(9) 



P = ^RcJ dyj{y)E{yr 

= ^Re Jdy Jdy'E{yra{y,y')E{y') 



(10) 



We consider the case that the length of the nanotube 
is larger than the mean free path. In this case, we can 
neglect effects of edges on the conductivity of the carbon 
nanotube and replace the conductivity with 



<j{y,y')^eo{y)eo{y'My-y'), 



(11) 



where cr(y—y') is the conductivity in infinitely long nan- 
otubes and Oo{y) is unity well inside the nanotube and 
should decay rapidly outside. The simplest choice is 



Z\2 



where 9(t) is the step function defined by 

m 



1 (i>0); 
{t<0). 



(12) 



(13) 



Actually, a nanotube has circumference L = 2t:R and the 
field should be smoothed out over the distance of the 
order of L. Therefore, we shall replace 0o(2/) with 



1 — erf 



(14) 



with effective edge width ^ = where the error function 
is defined by 



erf(i) 



-^'dt. 



(15) 




y=+//2 



y=-//2 



FIG. 1: A schematic illustration of a finite-length nanotube 
with length I. In actual numerical calculations, we shall con- 
sider a periodic array with period d [d^l). 



When the mean free path is much larger than the tube 
length, discrete energy levels are formed because of quan- 
tization of the electron motion and deviations from this 
approximation may appear. 

For simplicity, we consider the case that the Fermi level 
lies in the linear band 



e = ±^\k\ — ±hvp\k\, 



(16) 



where k is the electron wave- vector, v-p — is the 
Fermi velocity, and 7 = '\/3a7o/2 with the lattice con- 
stant a and the nearest-neighbor hopping integral 70 in 
graphene. Further, we shall employ a relaxation-time 
approximation to calculate the Boltzmann conductivity, 
which becomes an exact solution of the transport equa- 
tion in the single-channel case. As shown in Appendix 
A, we have the conductivity 

a(g,t^) =5vgs;^(-ic^^)[(-it^+^)"'+(wFg)^] \ (17) 

where and are the degeneracy for the valley and 
spin, respectively. The impurity and/or phonon scat- 
tering are characterized by relaxation time r, related to 
mean free path A given by 



K — V-pT. 



B. Periodic Array 



(18) 



For actual numerical calculations, we consider the ar- 
ray with period d as illustrated in Fig. [T] and seek the 
solution periodic with d. The electric field is expanded 
into the Fourier series 



i?(y)=^i?(G)exp(iGy), 



(19) 



with reciprocal lattice G — ^irj /d [j — 0, ±1, • • • ). Then 
the equation for the electric field becomes 



V Ug.g' - --,K'{G)ct{G,G')\e{G') = E,^,{G), 

^ — ' I ILO d J 

G' 



(20) 
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where 



Approximately, we have 



a(G,G') = / / ga(g,u;)^o(G-g)^o(G'-g)*, (21) 



and 



9o{q) - y 



(22) 



In terms of the Fourier coefficients, the power absorption 
of each nanotube is written as 



(23) 



G G' 



Because K'{q) increases in proportion to q for large g, 
the convergence of the solution of Eq. ((20)) requires very 
large values of G. This slow convergence corresponds to 
the singular distribution of the electric field associated 
with point-like polarization charges building up in the 
vicinity of the end points of the nanotube. Actually, the 
use of the nonlocal conductivity is quite effective in sup- 
pressing this singular behavior, but is not essential for 
determining resonance behavior particularly in the low 
frequency region. 



C. Single-Mode Approximation 

Let us consider the case in a uniform external field, 
i.e., £'ext(y) = £'oxt, which is valid for light with wave- 
length much larger than the tube length and polarized 
in the y direction. Polarization charge is induced at the 
ends of the finite length nanotube, as mentioned above. 
This charge is regarded almost as a point charge and the 
associated electric field decreases in proportion to l/r^, 
where r is the distance from the end point. In sufficiently 
long nanotubes, therefore, the field due to this induced 
charge is negligible in the most region of the nanotube 
where the absorption mainly takes place. If we neglect 
the presence of the ends, a single plasmon mode is in- 
duced along the nanotube determined by the wave num- 
ber corresponding to the frequency of the external field. 
For such a long-wavelength mode, the response is essen- 
tially determined by local conductivity criuj) = a{0,Lu), 
i.e., (j{y — y') — a{uj)d{'y — y'). Further, the main effect 
caused by the presence of the edges on this mode is to 
impose the boundary condition that the induced current 
should vanish at the edges, leading to the vanishing elec- 
tric field. 

Let Qi^ be the wave number of such a mode domi- 
nantly excited. Then, the electric field associated with 
this "bulk" mode is approximately given by 



E{y) « £;o[cos(Q^2/)-cos(Q„//2)], 



(24) 



with appropriate coefficient Eq. The second term in the 
bracket arises from the boundary condition, i?(±//2) = 0. 



I dy' J dy" K'{y-y')a{y', y")Eo [cosiQ^y") - cos{QJ/2)] 
« K'{Q^)a{u;)EoCOs{Q^y). (25) 
Thus, Eq. ([7]) gives two equations: 

"^"^--Q.) = 0, 

(26) 
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- Eqcos{QJ/2) = E^^t. 
The second equation gives 

E{y)^E.Jl 



cos(Q^//2). 
The first equation gives the mode frequency: 



1 



w[w+(i/t)] 



0, 



with 



(27) 



(28) 



(29) 



In the limit of an ideal nanotube r— > cx3, we have = 
zLQ with Q>0, satisfying 

a;{Qf^9^^2 (30) 
where vq is the approximate velocity of the mode 



4vp 



4 27rK7 



MQR)IoiQR). (31) 



Except in extremely short nanotubes with length not 
so much different from circumference L, the velocity is 
nearly constant because of logarithmically weak depen- 
dence on Q. We have e^/2Trfi-f < 0.2 and Kq{QR)Iq{QR) 
> 1, and therefore vq > v-p. In clean wires satisfying 
LUT^l, we have Q^i — Q + iQ' , with 



Vq 2vqT 2vq a 

In dirty wires lut<^1, on the other hand, we have 



Qu 



1 



V2vQ 



(l+i)- 



(32) 



(33) 



Figure [2] shows some examples of for tubes with 
1/ R = 2ttI/ L = 15Q0 (corresponding to a tube with length 
I w 1 fim and diameter 2R « 1.36 nm of the so-called 
(10,10) nanotube). We have used 70 = 3.0 eV and ko — 2.5 
corresponding to the bulk graphite. 42i In this example, 
we have roughly vq » Av-p. The velocity is larger than 
the electron velocity in the graphene or in metallic car- 
bon nanotubes, but is still much smaller than the light 
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In fact, in the vicinity of (5„, we have 
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FIG. 2: Examples of calculated plasmon dispersion in a tube 
with radius 7?/^ = 1/1500 for A// = 0.02, 0.1, and 0.5, where 
A is the mean free path. Solid and dotted lines show the 
real and imaginary parts of Q^, respectively. The real part is 
almost linear with the slope of vq/vf ~ 4 in the clean tube 
A// = 0.5. 



E{y) « £;ext(-l)"(cos 

■{Q-Qn)l 



1 , ' 

COS — — cosh ■ 



on 

1 tanh ^— 



(36) 



showing that the imaginary part exhibits a resonant in- 
crease at Qn following a Lorentzian form and the real 
part a Lorentzian multiplied by Q — Qn- Therefore, the 
absorption also exhibits a resonance proportional to 



( Q'l 

(cos — 



1 + [ujt) 



tanh 



(37) 



with (To = "■(OjO) = gvQsi^"^ l''^K)K and total number N 
of isolated nanotubes. 

This resonance behavior decreases with the increase of 
disorder or the decrease of mean free path A. The clear 
resonance disappears when cosh{Q'l/2) « sinh((3'//2), 
i.e., Q' ~ Ti/l. In the present example with vq ~ 4t!F, 
this condition becomes K/l ~ 0.04. In the dirty limit, 
the field remains small for small w or Q (» Q') until 
condition Q ^ Q' ^ tt/1 is satisfied and then becomes 
essentially the same as the external field except in the 
vicinity of the edges. As a result, the absorption takes 
a broad maximum around the corresponding frequency 
and decays with oj following the real part of the dynam- 
ical conductivity cr(w). 



III. NUMERICAL RESULTS 



velocity c because ~ c/300. This slow plasmon ve- 
locity justifies the use of Poisson's equation instead of 
full Maxwell's equation used for normal metal wires. 
The above results become essentially the same as that 
obtained in ref. |2^ based on the so-called Leontovich- 
Levin equation when we take the limit c— >oo. Plasmon 
modes have been theoretically studied in infinite one- 
dimensional organic conductors*'* and in semiconductor 
quantum wires 4^ 

In clean tubes satisfying condition cosh(Q'^/2) ^ 
smh{Q'l/2) and therefore cosh(Q'2//2) > | sinh((5'y/2)|, 
we have approximately 



E{v) 



, Qy , Q'y 

1 — cos — cosh 

2 2 



Ql ^Q'l . Ql . , Q'l\-^' 

X I cos — cosh ism — smh — 

2 2 2 2 



(34) 



except in narrow regions of y for which cos{Qy /2) w 0. 
This exhibits a resonance behavior at Q = Qn with 



g„=(2n + l)y (n = 0,l,---). 



(35) 



Figure [3] shows some examples of calculated absorption 
intensity in the nanotube the same as shown in Fig. [21 
The mean free path is chosen to be A/l = 0.5 (clean), 
0.1 (intermediate), and 0.02 (dirty) in these examples. 
The solid lines show numerical results and the dotted 
lines those obtained in the single-mode approximation 
discussed in the previous section. The positions corre- 
sponding to Q = Qn are denoted by thin vertical lines. 

The main resonance with n = occurs almost exactly 
at Q — Qo- For the higher- frequency resonance n=l, the 
peak position is slightly shifted to the higher frequency 
side. This shift becomes larger with the increase of n. 
In the single-mode approximation, on the other hand, 
the resonance occurs exactly at Qn- This shift is due to 
the appearance of strongly localized electric field associ- 
ated with induced charges at both ends of the nanotube. 
Apart from this peak shift, the single- mode approxima- 
tion works surprisingly well. 

Figure m shows the electric field at the center, y — 0, 
of the clean tube {A/l = 0.5) as a function of the fre- 
quency. The real and imaginary parts exhibit the res- 
onance behavior same as that discussed in the previous 
section within the single- mode approximation. In fact. 
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FIG. 3: The absorption in a nanotube with radius R/l — 
1/1500 as a function of frequency. A/l — 0.5, 0.1, and 0.02. 
The results in the single-mode approximation are shown by 
dotted lines. Corresponding wave number Qui in an infinitely 
long tube with A ^ oo is shown in the upper axis. Thin 
vertical dashed lines indicates wave number Qn — iv(2n+l) /I. 



the single-mode approximation has originally been ex- 
pected to be sufficient for the behavior of the field in the 
central part of the nanotube. 

Figure shows calculated electric- field distribution in 
a clean tube with A/l = 0.5. The dotted lines show the re- 
sults in the single-mode approximation. The most note- 
worthy feature is large electric field localized at both ends 
of the nanotube. This is due to the significant accumu- 
lation of induced charges at both ends of the tube. This 
field decays rapidly away from the tube edge and ap- 
proaches external field E^xt outside the tube. Well inside 
the tube, the field becomes approximately sinusoidal with 
wave vector Q of the plasmon mode determined by the 
frequency. In the low-frequency region corresponding to 
Q smaller than Qq, the electric field is screened out in a 
carbon nanotube by charges accumulated at both ends. 
At resonances, both real and imaginary parts of the field 
are considerably enhanced. 

Difference between the exact numerical results and 
those of the single-mode approximation becomes appar- 



FIG. 4: The real and imaginary parts of the electric field at 
the center of the nanotube as a function of frequency. The 
dotted lines represent those of the single-mode approximation. 



ent in the vicinity of resonance n = 2, i.e., Ql/'K = b. This 
arises mainly due to the shift in the peak frequency in 
the numerical result, not present in the single- mode ap- 
proximation. For larger n, the field distribution is more 
strongly affected by accumulated edge charges and extra 
restoring force tends to enhance the resonance frequency 
from that in the single-mode approximation. 

Figure[6]shows calculated electric-field distribution in a 
tube with A/Z = 0.1, for which Q'l^l. In this case broad 
resonance is recognized only in the imaginary part of the 
electric field and deviation from the single-mode results is 
less prominent. In a dirty tube with A/Z = 0.02 shown in 
Fig. [71 the field distribution is almost independent of the 
frequency except in the low frequency region where Q ~ 
Q' ^1^/1, as has been discussed in the previous section 
for the results in the single-mode approximation. 

Figure [8] shows the details of field distribution (imagi- 
nary part of E{y)) around the right end of the nanotube 
for (a) Ql/'K= \ and (b) 5. Polarized charges appear at 
the edge roughly in proportion to the conductivity and 
therefore their amount decreases with the decrease of the 
mean free path. A large electric field appears associated 
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FIG. 5: Calculated electric field distribution for varying Q in a clean wire with mean free path K/l — 0.5. The dotted lines 
represent the results of the single-mode approximation. Lines for increasing Ql/ir are shifted in the vertical direction, for each 
of which the horizontal axis at i? = is shown by a dashed line. 
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FIG. 6; Calculated electric field distribution in a wire with A/l — 0.1. 
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FIG. 7: Calculated electric field distribution in a dirty wire with A./1 = 0.02. 



with these charges and decays slowly in proportion to 
{y—\l)~'^ outside of the nanotube. This outside field 
does not contribute to the power absorption, however. 

On the other hand, the electric field due to edge 
charges rapidly decays inside the nanotube due to strong 
screening effect and becomes negligible in the most region 
of the nanotube where the absorption mainly takes place. 
In fact, only small difference appears in the absorption 
power given by the solid and dotted lines in Fig. [31 This 
shows that the strong screening of the field due to polar- 
ized edge charges is the main reason that the single-mode 
approximation based on a plasmon mode in an infinitely 
long nanotube has turned out to work surprisingly well 
in spite of the presence of edges. 



IV. DISCUSSION AND CONCLUSION 

Carbon nanotubes used in experiments'^*^^ are usually 
several micrometer long and therefore longer than the 
mean free path limited by impurity and/or phonon scat- 
tering at room temperature. Recent progress in exper- 
iments has achieved short nanotubes with sub-micron 
size^i^ili^^^ In such short nanotubes, discrete energy 
levels are formed and the level spacing can exceed the 
broadening due to disorder. Then, we should seriously 
consider effects of discrete energy levels to discuss opti- 
cal absorption of finite-length nanotubes. This problem 



is out of the scope of this paper and left for a future 
study. 

When a nanotube is tilted from the direction of the ex- 
ternal field, the field component parallel to the axis is ef- 
fective in the absorption and the resonance frequency re- 
mains unaffected. When many nanotubes are distributed 
at random within a plain, the absorption intensity is re- 
duced by factor 1/2 due to the average over directions. 
Recently, a monotonic dependence on the direction is re- 
ported in a mat of aligned nanotubes. For the field per- 
pendicular to the axis, there is no significant absorption 
in the relevant frequency range considered heTeM^^^ In 
fact, the dynamical screening or the depolarization effect 
due to induced charges should be consideredii^ and res- 
onance absorption appears at much higher energies when 
exciton effects are properly included in interband optical 
transitions.'^^!-^ 

In actual absorption experiments, a bundle of nan- 
otubes may be used. When a bundle contains N metallic 
nanotubes with same length, the Coulomb kernel given 
by Eq. ([5]) is effectively multiplied by N . As a result, 
the frequency and the velocity of excited plasmon mode 
are multiplied by ^/N . Semiconducting nanotubes can 
be contained in a bundle, but usually have very small 
conductivity and irrelevant, although it is known that 
carbon nanotubes tend to be naturally doped due to sur- 
rounding particles such as oxygen 

Usually, experiments are performed for a film where 
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FIG. 8: A close-up view of imaginary part of electric field 
distribution around the right end of nanotubes for (a) Ql/n = 
1 and (b) 5. 



the density of nanotubes is sufficiently small for the pur- 
pose of avoiding electrical contacts among them. In this 
case, electric field induced by neighboring nanotubes can 
be considered in a dipole approximation. The dipole 
field decays as with the increase in distance r. In 
film-like systems, this dipole field of distant nanotubes 
does not give rise to significant contribution in con- 
trast to three-dimensional systems. Effects of surround- 
ing nanotubes were previously considered by numerically 
solving Hallen's equation for infinite planar arrays of 
nanotubes i^i^l The results showed that interaction be- 
tween neighboring tubes causes broadening and shift of 
resonances only when the distance becomes of the order 
of the tube diameter. If an end of a nanotube touches 
another nanotube, the absorption may be significantly 
modified because of large charges accumulated at the 
end. Such a situation should be avoided in experiments. 

In general, a mini-gap opens around the Fermi energy 
even in so-called metallic nanotubes for various reasons. 
For example, the nonzero curvature causes a band gap 
depending on the chirality,^^ which can be understood in 
terms of an effective Aharonov-Bohm fiux within the k-p 
schemei^i^i^ A mechanical deformation causes band- 
gap modification;^!^ which can also be understood in 
terms of fluxi^i^ The amount of the gap depends on 
the chirality, radius, environment, etc., and is typically 
smaller than room temperature. There have been some 
reports suggesting the observation of broad absorption 
due to inter-minigap transitions . '^^^'^^ Because this mini- 
gap transition is in the same frequency region, detailed 
study on the length dependence is required for the pur- 
pose of identifying the absorption due to the finite- length 
origin considered here. 

In summary, we have calculated electric-field distribu- 
tion and absorption intensity of a finite-length nanotube 
in oscillating electric field. The results show that the 
main resonance corresponding to excitation of the funda- 
mental plasmon mode with wave vector Qo = 'k/1 is quite 
robust except in very dirty tubes. For higher-frequency 
resonances, the electric field associated with induced edge 
charges starts to be mixed and tends to shift resonances 
to higher frequencies. Overall resonance behaviors can 
be reasonably well described by the single-mode approx- 
imation in which effects of induced edge charges are com- 
pletely neglected. 
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APPENDIX A: NONLOCAL CONDUCTIVITY 



The induced current 



We start with a Boltzmann transport equation 

+ Vjk ' -eE{y,t) 



dfjk{y,t) dfjkiy,t) ldfjk{y,t) 



dt dy 'h dk 

= -^Wfk'jk[fjk{y,t)-fj'k'{y,t)], (Ai) 

j'k' 

where j denotes subbands, k the wave vector in the axis 
(y) direction, E{y, t) is the apphed electric field, Vjk = 
dejk/fidk is the velocity, and Wj'k'jk is the scattering 
probability between states jk and j'k'. We write the 
distribution function fjk as the sum of the equilibrium 
distribution function f{ejk) and the deviation due to the 
applied field gjk, i.e., fjk{y,t) = f{sjk)+gjk{y,t)- In the 
limit of a weak applied field, the above is approximated 
by 



j{q,uj) = gyQs ^{-e)vjkgjk, (AlO) 

jk 

can be expressed in terms of the conductivity a{q, w) as 

j{q,oj) = a{q,co)E, (All) 

with 



a{q,u)) = j (- 



de 



and 



a{q,uj,e)de, 



(A12) 



dgjk{y,t) dgjk{y,t) 
1- Vjk ^ h evj 



dt 



dy 



= - ^ Wfk'jk [gjk{y, t)-gyk' {y, t)]. 

j'k' 

Let 

E{y,t) = Eexp{-iu)t+iqy). 
Then, we can set 

9jk{y,t) = gjkexp{-iu}t+iqy), 
and rewrite the transport equation as 



dsjk 
(A2) 

(A3) 
(A4) 



a{q, u), s) = 5fv5s Y] \ - iu)+iqvjk+-r^] e'^v'^J{e-ejk). 

(A13) 

Adding the contributions of positive and negative values 
of A;, we can rewrite the above as 



a{q,uj,s) = g^gs^ e'^\v-jk\'^ 



jk>0 



1 \2 



-lU! 



r{e) 



T{e) 
6{s-Sjk). 



-q mk\ 



(A14) 



We shall consider the energy region within the linear 
band. Substituting the explicit expression Vk = VY=^/h 
and ek='yk, we have 



- ~d^) o-(g, LO, e) = g^g, 



2 



1 



1\2 72^2 



= - y~] Wjfk'jklgjk-gj'k']- 
j'k' 

The solution can be written as 

dfisjk) ' 



9jk 



(A5) 



(A6) 



with T = T{e) for simplicity. 

The poles of the conductivity are given by 

huj = ±79 — 1 — , 

T 



(A15) 



(A16) 



with 



[-iu>i-iqvjk]gjk+evjkE = - ^ Wfk' ,jk[9jk-9j'k']- (A7) 

j'k' 

This constitutes a set of linear equations determining gjk 
for each energy and therefore can be solved exactly. In 
the following, however, we shall use a simplest relaxation- 
time approximation, which become exact for a single 
channel case. 

The right hand side is approximated by 



which corresponds to the subband dispersion in the limit 
of T^oo, as is expected. In the limit of small 1/r and w, 
the denominator is a function of Dq^—ioj as is expected, 
where the diffusion constant is given by 



(A17) 



^Wj'k',jk\gjk{y,t)-gj'k'{y,t)] = --7^gjk{y,t). 

T\^jk) 

(A8) 



In the high-frequency limit ujt^I, on the other hand, 
the cut-off wave vector q^, where the conductivity devi- 
ates from the local conductivity a{Lo) = a{Q,uj), becomes 
Uf9c~w or 



j'k' 



Then, we immediately have 



2-K \ L J ' 



(A18) 



9jk 



- iLij+iqVjk+' 



1 



T{ejk) 



evjkE 



dsjk 



■ (A9) 



which should be much smaller than unity under the usual 
conditions. 
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